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Abstract 

We analyze the Bell polynomials B n {x) asymptotically as n —* oo. We obtain asymptotic 
approximations from the differential-difference equation which they satisfy, using a discrete 
version of the ray method. We give some examples showing the accuracy of our formulas. 
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1 Introduction 

The Bell polynomials B n (x) are defined by [T] 

n 

B n {x) = Y,S n k x\ 71 = 0,1,..., 

where S% is a Stirling number of the second kind [2J 24,1,4]. They have the generating function 

°° t n 

J> n (*)- = ex P [x(e'-l)], (1) 



n=0 
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from which it follows that 

B (x) = 1 (2) 

and 

B n+1 (x) = x[B' n (x)+B n (x)], n = 0,l,.... (3) 

The asymptotic behavior of B n (x) was studied by Elbert [3], [lj and Zhao [5], using the saddle 
point method and ([1]). In this paper we will use a different approach and analyze Q instead of ([1]). 
The advantage of our method is that no knowledge of a generating function is required and therefore 
it can be applied to other sequences of polynomials satisfying differential-difference equations [6], 

2 Asymptotic analysis 

To analyze ([3]) asymptotically as n — > oo, we use a discrete version of the ray method [Hj. Replacing 
the anszat 

B n (x) =e- n F(ex,en) (4) 

in ([3]), we get 

(dF \ 

with 

u = ex, v = en (6) 
and £ is a small parameter. We consider asymptotic solutions for ([5]) of the form 

F(u, v) ~ exp [e~V ( M ) v )] K{u, v ), (7) 

as e — > 0. Using (|7j) in (JSJ) we obtain, to leading order, the eikonal equation 

e q -u(p + l) = (8) 

and the transport equation 

dK l^± R _ uex (_^f)^L_ Q (9) 
dv 2 dv 2 U 6 ^ \ dv J du 

where 

dx' dv' 

The initial condition (j2J), implies 

V>(u,0) = 0, J fsT(M,0) = l. (11) 

To solve © we use the method of characteristics, which we briefly review. Given the first order 
partial differential equation 

$ (u,v,ip,p,q) = 0, 
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with p, q defined in (fTUI) . we search for a solution ip{u,v) by solving the system of 
equations" 

du 9$ dv d$ 

dt dp ' dt dq' 

dp d$ d$ . dq d$ d$ 
^ dt du ^ dip' ^ dt dv dip 1 

dt dp dq ' 

where we now consider {u, v, ip,p, q} to all be functions of the new variables t and s. 
For ([8]), we have 

$ (u, v, ip,p,q) = e q + p - 2u 
and therefore the characteristic equations are 

u + u = 0, v = e q , p — p=l, q = 

Solving ffl2l) . subject to the initial conditions 

u{0,s) = s, v{0,s) = 0, p(0,s) = B(s) - 1, 

we obtain 

u = se~\ v = Bst, p = Be 1 - 1, q = In (Bs) 

where we have used 

= d\ t=0 = e 9(0 ' s) - sB. 

From ( TTTT) and ( fl3l) we have 

^(0,s) = 0, tf(0,s) = l, 

which implies 

= (0, s) = p(0, s)^« (0, s) + g(0, s)^-v (0, s) 
= (5 - 1) x 1 + In (5s) x = B - 1. 

Thus, 

U = se~*, v = st, p — e — 1, g = In (s) . 
The characteristic equation for ^ is 

ip = pii + qv = (e* — l) (— se~') + In (s) s, 

which together with (TH1) gives 

V> (t, s) = s (1 - t - e _t ) + In (s) st. 

We shall now solve the transport equation (jUJ). From f|T5|) . we get 

9t £e* 9t 1 9s e* 9s 1 

9m = ~s(t+ 1)' 9v = s(t + 1)' 9m = t + V dv = 7+1 



and therefore, 



d 2 ijj dq dq dt dq ds 1 



dv 2 dv dt dv ds dv s (t + 1) 
Using (fT71) - (|T8l) to rewrite ([9]) in terms of t and s, we have 

1 



K 



2(t + l 



-K = 



with solution 



where we have used (fl4"|) . 

Solving for i, s in (TT5I) . we obtain 



K{t,s) 



Vt + T' 



V 



t = UN - , s = 

LW (s) 



where LW (•) denotes the Lambert-W function [9], defined by 

LW (z) exp [LW {£)} = z. 
Replacing (j2"0|) in (jlBjl and (IT91. we get 



K(u, v) 



+ v\n 



i 



LW(f)+l 



and from (J7|) we find that 

t>) ~ exp 



vie v , 
+ - In 



LW (~) 



lw(3 



,u + v N 



Lw©+i 



as e — > 0. Using (jSJ) and (j2ll) in (j4j), we conclude that 



B n (x) ~ exp 



LWg) 



+ n In 



n 



LW(*) 



Lw(a+i 



as n-> oo. 



(18) 



(19) 



(20) 



(21) 



(22) 



Remark 1 The function LW (z) /ias two real-valued branches for —e 1 < z < 0, denoted by LWo (z) 
(£/je principal branch of UN) and LW_i (z) , satisfying 

LWo : [-e"\ 0) -+ [-1, 0) , LW_! : [-e~\ 0) -> (-oo, -1] , 
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with 

LW (-e" 1 ) = -1 = LW_! (-e" 1 ) . 
For z > 0, LW (z) has only one real-valued branch 

LW : [0, oo) -> [0, oo) 

and for z < — e -1 , LWo (z) and LW_i (z) are complex conjugates. Therefore, for [2^1 to be well 
defined, we need to consider three separate regions: 

1. An exponential region for x > or x < —en. Here we have 

B n (x) ~$ n (x;0), n^oo, (23) 

where 



$ n (x; k) = exp t /\ , + rcln 



(as + n) 



lw, o + 1 



yln oscillatory region for —en < x < 0. In this interval, 

B n (x) ~ $„(as;0) + $ n (x;-l), n -> oo. (24) 

In FigureU\ (a) we plot B 5 (x) and the asymptotic approximations If23\) (+++) and (24\ ) (ooo), 
all multiplied by e~' x ' for scaling purposes, in the interval (—10, 10) . We see that our formulas 
are quite accurate even for small values of n and that the transition between (23}) and [Efy is 
smooth. 

3. A transition region for x ~ —en. We will analyze this region in the next section. 

In Figured] (b) we plot B 5 (x) and fl23l (+++) and fl24|) (ooo), all multiplied by e x , in the interval 
—20, 0) . We observe that the approximations ( 1231) and ( 1241) break down in the neighborhood of 
-eh ~ -13,59. 



2.1 The transition region 

When x = —en, the quantity LW (-) +1 vanishes and (1231 is no longer valid. To find an asymptotic 
approximation in a neighborhood of —en, we introduce the stretched variable (3 defined by 

x = - en -(3n^, (3 = 0(1). (25) 

For values of z close to zq = — e" 1 , the Lambert-W function can be approximated by [9j (4.22)] 

LW (z) ~ -1 + y/2e (z-zq) - ~e (z - z ) + ^2e 3 (z - z ) 3 , z -e" 1 . (26) 
Using (l2~5]) in (|26|) . we have, 



(a) (b) 
Figure 1: A comparison of the exact (solid curve) and asymptotic (ooo), (+++) values of B${x). 



n 



+ nln 



n 



{x + n) } ~</?(M, P^O, 



Hence, 

for k = 0,1 with x = —en — /3ns and 

p (/?, n ) = (-l) n exp | [ln(n) + e - 2] n - (e" 1 - l) /3n^ } . 
We now consider solutions for ([3]) of the form 

B n {x) = <p(j3,n)A(j3) = <p 

for some function A (/?) . Replacing (j2"9"|) in (j3J) and using (T25]) we obtain, to leading order 

A" - 2e~ 3 /3A = 0, 

with solution 





X\ 2 


A 




[-(eH 


— ) n 3 ,n 


- ( e+ D nl . 




nJ 







A (/?) = dAi (2 V 1 /?) + C 2 Bi ^e" 1 /? 



(28) 



(29) 



(30) 



where Ai (•) , Bi (•) are the Airy functions. 

To determine the constants C x , C 2 in flUJ, we shall match (J23J with (J2U). Using (j25) and (T27D 
in fl23l). we have 



5 n (x) 



<^(/3,n)exp ^— |V2e _ t/5iJ (2e _1 /3)* 



(31) 



On the other hand, the Airy functions have the well known asymptotic approximations [2J (10.4.59, 
10.4.63)] 

Ai(z) ~ 2^ 6XP (~3^) Z ~^ °°' 

1 /2 s\ _i 
Bi (z) ~ — p= exp -z 2 z 4 , 2 — > oo 
/7T V 3 / 



and therefore we conclude that 



Ci = 0F2fnS, C 2 = 0. (32) 
Replacing fl30|) and fl32|) in fl29l) . we find that for x ~ —en, we have 

B n (x) ~ v^^rii^ (/3, n) Ai ^e -1 /?^ , n — > oo. 

This concludes the asymptotic analysis of B n (x) for large n. 
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